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Abstract 

We list in increasing order - - {± ± , §, §, ±§£, ^, |, f, ff, 2 - |§f , ±| , l} - - a number 
of exact two-qubit Hilbert- Schmidt (HS) separability probabilities, we are able to compute. 
Each probability corresponds to a specific scenario — a class of 4 x 4 density matrices (p) 
with an m-subset (m < 6) of its six off-diagonal pairs of symmetrically located entries set to 
zero. Initially, we consider only scenarios in which the (6 — m) non-nullified pairs are real, but 
then permit them to be complex (as well as quaternionic) in nature. The general analytical 
strategy we implement is based on the Bloore density matrix parameterization (J. Phys. A, 9, 
2059 [1976]), allowing us to conveniently reduce the dimensionalities of required integrations. 
For each scenario, we identify a certain univariate "separability function" S sce nario( l ')i where 
v = pWpzl ■ The integral over v 6 [0, oo] of the product of this function with a scenario-specific 
(marginal) jacobian function J S cenario{v) yields the HS separable volume (V^f). The ratio of 
Vsep to the HS total (entangled and non-entangled) volume gives us the HS scenario-specific 
separability probability. Among the possible forms that we have so far determined for 5 scenar j (^) 
are piecewise combinations of c, C\J~u and cu for v G [0, 1] and (in a dual manner) c, ^= and ^ 
for v £ [1, oo]. We also obtain bivariate separability functions S®*® nario {vi,V2) in the <\ub\t- qutrit 



case, involving 6x6 density matrices, having ratio variables, v\ = pllp55 and uo = p22p66 . 

' o ' 1 P22P44 Z P33P55 

Additionally, we investigate parallel two-qutrit and three-qubit problems. Further still, we find 
some analytic evidence for the relevance of beta functions to the two-qubit separability function 
problem, while we previously (Phys. Rev. A [2007]) had found numerical evidence of such a nature. 
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I. INTRODUCTION 

In the abstract to their recent comprehensive review, the Horodecki family note that while 
quantum entanglement is "usually fragile to environment, it is robust against conceptual and 



4 



mathematical tools, the task of which is to decipher its rich structure" In the study 
below, we certainly encounter such robustness and pursue the task of surmounting it — with 
some successes. 

Given a generic class composed of composite quantum systems, the question of the "rel- 
ative proportion" of entangled and non-entangled states in that class was apparently first 
raised by Zyczkowski, Horodecki, Sanpera and Lewenstein (ZHSL) in a much-cited paper 
Q]. They gave "three main reasons" — "philosophical", "practical" and "physical" — upon 
which they elaborated, for pursuing the topic. The present author, motivated by the ZHSL 
paper, has investigated this issue in a number of settings, using various (monotone and non- 
monotone) measures on quantum states, and a variety of numerical and analytical methods 



m (cf. 



11 



12, 



141] ). Though the problems are challenging (high- 



dimensional) in nature, many of the results obtained in answer to the ZHSL question in 
these various contexts have been strikingly simple and elegant (and/or conjecturally so). 
The new results below certainly fit into such an interesting, appealing pattern. 

Specifically here, we develop an approach recently presented in 15]. This was found to 
be relatively effective in studying the question posed by ZHSL, in the context of two-qubit 
systems (the smallest possible example exhibiting entanglement), endowed with the (non- 
monotone [if]]) Hilbert- Schmidt (HS) measure 17], inducing the flat, Euclidean geometry 
on the space of 4 x 4 density matrices. This approach is| exploits two distinct features of 



a form of density matrix parameterization first discussed by Bloore [18l |. These properties 
allow us to deal with lower- dimensional integrations (more amenable to computation) than 
would otherwise be possible. We further find that the interesting advantages of the Bloore 
parameterization do, in fact, carry over — in a somewhat modified fashion — to the qubit- 
qutrit (sec. IIIip . qutrit-qutrit (sec. IIVI) and qubit-qubit-qubit (sees. IVl and IVIj) domains. 
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A. Bloore (off-diagonal-scaling) parameterization 



We shall, first, consider the 9-dimensional convex set of (two-qubit) 4x4 density matrices 



with real entries, and parameterize them — following Bloore 



18) (cf. [li, p. 235]) — as 



Pn 



Z 12\J ' PllP22 Zis^/pu P33 ^14a/PhP44 ^ 



Z\2^ Pllp22 P22 ^23a/P22P33 z 2i a//?22P44 

z 13\> r Pllp33 Z 23\J 'P22P33 P33 2; 34 V / P33P44 



(1) 



/ 



\ ^14yPllP44 ^24a/P22P44 ^34\/P33p44 P44 

One, of course, has the standard requirements that > and (the unit trace condition) 
TtiPu = 1. Now, three additional necessary conditions (which can be expressed without using 
the diagonal entries, due to the pa > stipulation) that must be fulfilled for p to be a density 
matrix (with all eigenvalues non-negative) are: (1) the non-negativity of the determinant 
(the principal 4x4 minor), 



(^34 — Z 12 + ^ (Zu (^24 — z 23 z 3i) + %13 ( z 23 ~ Z 2A Z 34,)) z 12 



'23 



'21 



^34 + 



(2) 



+ ^14 ( Z 23 — l) + Z 13 ( Z 24 ~~ l) + 2223^24^34 + 2>Zl3 z U (^34 — 223^24) + 1 > 0; 

(2): the non-negativity of the leading principal 3x3 minor, 

— z\ 2 + 2z 13 2:232l2 — z i3 — ^23 + 1 > 0; 



(3) 



and (3): the non-negativity of the principal 2x2 minors (although actually only the i = j = 1 
case is needed, it is natural to impose them all), 



(4) 



As noted, the diagonal entries of p do not enter into any of these constraints — which taken 
together are sufficient to guarantee the nonnegativity of p itself — as they can be shown to 
contribute only (cancellable) non-negative factors to the determinant and principal minors. 
This cancellation property is certainly a principal virtue of the Bloore parameterization, 
allowing one to proceed analytically in lower dimensions than one might initially surmise. 
(Let us note that, utilizing this parameterization, we have been able to establish a recent 
conjecture of Mansson, Porta Mana and Bjork regarding Bayesian state assignment for three- 
evel quantum systems, and, in fact, verify our own four-level analogue of their conjecture 



20, eq. (52)] |2l|.) 



Additionally, implementing the Peres-Horodecki condition [221 . l23l . |24[ requiring the non- 
negativity of the partial transposition of p, we have the necessary and sufficient condition 
for the separability (non-entanglement) of p that (4): 

v {z\± ~ l) A2 + (uz 13 z u + z 23 z 2i - \fv (z u z 23 + z 13 z 24 ) 2:34) z l2 - z\ 3 - vz 2 3A + u+ (5) 

+V ((2:34 - l) zf 3 - 2z U Z 23 Z 24 Z 13 ~ ^24 + Z U { Z 23 ~ U )) + \z 13 Z 23 + VZ U Z 24 ) Z U > 0, 

where 

2 Pll^44 / c \ 

v = p = , (6) 

P22P33 

being the only information needed, at this stage, concerning the diagonal entries of p. (It 
is interesting to contrast the role of our variable u, as it pertains to the determination 
of entanglement, with the rather different roles played by the concurrence and negativity 



25|, [26].) We have vacillated between the use of v and p as our principal variable in our two 
previous studies [la, l27|. In sec. IVHl we will revert to the use of p, as it appears that its use 
avoids the appearances of square roots which definitely seem to impede certain Mathematica 
operations. 

So, the Bloore parameterization is evidently even further convenient here, in reducing 
the apparent dimensionality of the separable volume problem. That is, we now have to 
consider, initially at least, only the separability variable v rather than three independent 
(variable) diagonal entries. (This supplementary feature had not been commented upon by 
Bloore, as he discussed only 2x2 and 3x3 density matrices, and also, obviously, since the 
Peres-Horodecki separability condition had not yet been formulated.) The (two variable - 
Vi,v 2 ) analogue of (jSJ) in the 6x6 (qubit- qutrit) case will be discussed and implemented 
in sec. IHII Additionally still, we find a /bwr-variable counterpart in the 9x9 qutrit-qutrit 
instance [sec. [TV], and three- variable counterparts in two sets of qubit-qubit-qubit analyses 
(sees. IVl and IVll) . (The question of whether any or all of these several ratio variables are 
themselves observables would appear to be of some interest.) 

eqs. (3)-(5)], we expressed the conditions (found through application of the 



In 



15 



cylindrical algebraic decomposition" 28]) that - - in terms of the arbitrary 



9-dimensional 4x4 real density matrix p must fulfill. These took the form, 

Zl2, z 13, Zu G [— 1, 1], 023 G [Z 23 , Z 23 ],zu G [Z 24 , Z 24 \,Zu G [Z u , Z34], (7) 
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where 

Zf z = z 12 z 13 ± \Jl — z\ 2 \j\ — zf 3 , Z^ = zi 2 z u ± yjl — z\ 2 \J 1 — zf 4 , (8) 

y± _ ^13^14 — Z12Z14Z23 — Z 12 Z 13 Z 24: + ^23^24 i s 
34 ~~ 1 _ ^2 ) 

and 



— 1 + 2^ 2 + zf 3 — 2z 12 2l3^23 + Z 23\/ ^ + 2 12 + ^14 ~~ ^ z i 2 Zi4Z 2 4 + £24. (9) 



In his noteworthy paper, Bloore also presented [181, sees. 6,7] a quite interesting discussion 
of the "spheroidal" geometry induced by his parameterization. This strongly suggests that it 
might prove useful to reparameterize the Zij variables in terms of spheroidal-type coordinates. 
Closely following the argument of Bloore — that is, performing rotations of the (^13,^23) 
and (z\4, z 2i ) vectors by j and recogizing that each pair of so-transformed variables lay in 



ellipses with axes of length ^Jl ± z\ 2 — we were able to substantially simply the forms of 
these conditions (©-(ED). 

Using the set of transformations (having a jaoobian equal to 



z 13 — ► 7= , (.1UJ 



^23 



z 14 



Z24 



(Vi- 


- Z12 COS (#1) 


+ sin {fix) 


sjzvi + 1 ) V7172 + 1 


(sin ( 








Oi ) y/z 12 + 1 


— COS (#1) 


V 1 - Z X 2 ) V1112 + 1 






V2 




(VT 


- Zi 2 COS (9 2 ) 


+ sin (9 2 ] 


V^12 + l) \/7l + 72 


(sin 






02 ) VZ12 + 1 


- cos (9 2/ 


V 71 - z i 2 ) V7i + 72 



V2 



cos (6>i - 2 ) V71 + 72V7172 + 1 

^34 - ► ^34 



'72 

one is able to replace the conditions (©-(El)) that the real two-qubit density matrix p - 
given by ([1]) — must fulfill by 

71 e [0,1]; 7 2 e[7i,-]; ^34 e [-71,71]; 2 12 e [-M]; 0i,0 2 e[o,24 (11) 

7i 

We became aware of this set of transformations at a rather late stage of the research reported 
here, and have not been able so far — somewhat disappointingly — to exploit it in regards 
to the HS separability-probability question. So, the results reported below rely essentially 
upon the conditions (((7j)-(jSJ)) and the parameterization in terms of the z^s of Bloore.. 
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B. Previous analysis 



In 



151 ]. we studied the four nonnegativity conditions (as well as their counterparts 



having completely parallel cancellation and univariate function properties - - in the 15- 
dimensional case of 4 x 4 density matrices with, in general, complex entries) using numerical 
(primaril y q uasi-Monte Carlo integration) methods. We found a remarkably close fit to the 
function [15;, Figs. 3, 4], 



entering into our formula, 



4+ 5^)KHMH' 



(12) 



V. 



HS 

ep/real 



2 / Jreal{v)S rea i{v)dv 



Jreal(v)S rea ,l(v)dv, 



(13) 



or the 9-dimensional Hilbert-Schmidt separable volume of the real 4x4 density matrices 
9)1 . Here, B denotes the (complete) beta function, and B v the incomplete beta 



151 . eq. 



function 



29), 



B v (a,b) = / w^il-wf^dw. 

I — I 

Additionally [lj, eq. (10)], 

v 3 ' 2 (12 {y{y + 2) {y 2 + 14i/ + 8) + 1) log (y/u) - 5 (5z/ + 32z/ 3 - 32// - 5)) 



(14) 



(15) 



15 



3780(z/-l) 9 
Fig. 1]) jacobian function resulting from the 



is the (highly oscillatory near v 
transformation to the v variable of the Bloore jacobian (Uf =1 pa)^. (Perhaps we should 
refer to Jreaiiy) as a marginal jacobian, since it is the result of the integration of a three- 
dimensional jacobian function over two, say pn and P22, variables.) 



C. Computational limitations 

Although we were able to implement the three (six-variable) nonnegativity conditions 
(d2J), and (j4j)) exactly in Mathematica in [jjj], for density matrices of the form (CQ), 
we found that additionally incorporating the fourth Peres-Horodecki (separability) one (jHJ) 
- even holding v fixed at specific values — seemed to yield a computationally intractable 



problem. In fact, after the completion of 15j, we consulted with A. Strzebonski (the resident 



expert on these matters at Wolfram Inc.), and he wrote in regard to our problem that "It 
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looks like the [nine- dimensional four-condition separable real density matrix] problem is 
well out of range for CAD [the cylindrical algorithmic decomposition [3^]]. The algorithm 
is doubly exponential in the number of variables. Six variables is a lot for CAD, so only 
very, very simple systems with six variables can be solved. Adding one more inequality of 
total degree 4 makes a huge difference. After an hour the algorithm is still in the projection 
phase at 4 variables (it needs to go down to univariate polynomials) and the projection 
polynomials already are huge: the last resultant computed has degree 60 and 9520 terms, 
and the 3- variable projection set already has 890 such polynomials..." 



D. Research design and objectives 

In light of the apparent present computational intractability in obtaining exact results in 
the 9-dimensional real (and a fortiori 15-dimensional complex) two-qubit cases, we adjusted 
the research program pursued in We now sought to determine how far we would have 
to curtail the dimension (the number of free parameters) of the two-qubit systems in order 
to be able to obtain exact results using the same basic investigative framework. Such results 
- in addition to their own intrinsic interest — might help us understand those previously 
obtained (basically numerically) in the full 9-dimensional real and 15-dimensional complex 



cases 



15|. 



To pursue such a strategem, we chose to nullify various m-subsets of the six symmetrically- 
located off-diagonal pairs in the 9-parameter real density matrix (111) , and tried to exactly 
implement the so-reduced non- negativity conditions (([2]), (j3J), (j3J) and (jSJ)) — both the first 
three (to obtain HS total volumes) and then all four jointly (to obtain HS separable volumes). 
We leave the four diagonal entries themselves alone in all our analyses, so if we nullify m 
pairs of symmetically-located off-diagonal entries, we are left in a (9-m)-dimensional setting. 
We consider the various combinatorially distinct scenarios individually, though it would 
appear that we also could have grouped them into classes of scenarios equivalent under local 
operations, and simply analyzed a single representative member of each equivalence class. 

We will be examining a number of scenarios of various dimensionalities (that is, differing 
numbers of variables parameterizing p). In all of them, we will seek to find the univariate 
function S scenar i {v) (our primary computational and theoretical challenge) and the constant 
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^scenario-; SUCH that 



oo 







^sep / scenario ~ / ' scenarioiy) Cf scenarioiy) dl> , (16) 

and 

T ' 'tot/ 'scenario ~ ^scenario 

Given such a pair of volumes, one can immediately calculate the corresponding HS separa- 
bility probability, 

V HS , 



^tot / scenario ^scenario I 3 scenarioiy) dv '. (17) 



pHS _ se P/ scenario 

sep/scenario \/HS ' V / 



tot /scenario 



Let us note that in the 



ull 9-dimensional real and 15-dimensional complex two-qubit 
cases recently studied in [l5 ] , it was quite natural to expect that iS rea /(z/) = <S rea i(^) (and 
Scompiex{v) = S complex (£) ) • But, here, in our /ower-dimensional scenarios, the nullification of 
entries that we employ, breaks symmetry (duality), so we can not realistically expect such a 
reciprocity property to hold, in general. Consequently, we adopt the more general, broader 
formula in (fT3]) as our working formula ffl6l) . 



II. QUBIT-QUBIT ANALYSES 

To begin, let us make the simple observation that since the partial transposition operation 
on a 4 x 4 density matrix interchanges only the (1,4) and (2,3) entries (and the (4,1) and 
(3,2) entries), any scenario which does not involve at least one of these entries must only 
yield separable states. 



A. Five nullified pairs of off-diagonal entries — 6 scenarios 

1. 4- dimensional real case — = || 

There are, of course, six ways of nullifying five of the six off-diagonal pairs of entries of p. 
Of these, only two of the six yield any non-separable (entangled) states. In the four trivial 
(fully separable) scenarios, the lower- dimensional counterpart to S rea x(v) was of the form 

Sscenarioiy) ^scenario 2. 

In one of the two non-trivial scenarios, having the (2,3) and (3,2) pair of entries of p left 
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intact (not nullified), the separability function was 



2Jv < v < 1 
%,s)](") H ■ (19) 

2 i/ > 1 



(It is of interest to note that B v (\,l) = 2y/u, while in [15(, we had conjectured that S rea i(u) oc 
B v {\, y/3) and S complex {u) oc £„(^, ^§).) 

In the other non-trivial scenario, with the (1,4) and (4,1) pair being the one not nullified, 
the separability function was — in a dual manner (mapping f(y) for v G [0, 1] into /(-) for 
v G [1, oo]) — equal to 

{2 < v < 1 
(20) 

In both of these scenarios (having c scenar i = 2) for the total (separable and non-separable) 
HS volume, we obtained V t * s = ^ w 0.0654498 and V s f/ = ^ ~ .0385531. The cor- 
responding HS separability probability for the two non-trivial (dual) scenarios is, then, 
f£ w 0.589049. 

16 



5. 5- dimensional complex case — P^p = 3 

Although our study here was initially intended to concentrate only on 4 x 4 density 
matrices with solely real entries, at a later point in our analyses, we returned to (first) the 
m = 5 case, but now with the single non-nullified pair of symmetrically-located entries being 
complex in nature (so, obviously we have five variables/parameters in toto to consider, rather 
than four). 

Again, we have only the same two scenarios (of the six combinatorially possible) being 
separably non-trivial. Based on the (2,3) and (3,2) pair of entries, the relevant function 
(with the slight change of notation to indicate complex entries) was 

f-KV < V < 1 
(21) 
7T V > 1 



and, dually, 



7T < v < 1 

«%;4])]W = \ ■ (22) 

^ v > 1 

v 
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So, the function y/v, which appeared ( (TT9I) . (1201 ) in the corresponding scenarios restricted to 
real entries, is replaced by v itself in the complex counterpart. (We note that B u (l, 1) = v.) 

For both of these complex scenarios, we had V t ^ t s = and V^f = for a particularly 
simple HS separability probability of \. 



3. 7- dimensional quaternionic case — P^ s 



sep 



J_ 

10 



Here we allow the single pair of non-null off-diagonal entries to be quaternionic in nature 



31 



32| [33|, sec. IV]. We found 



and, dually, 



2v 2 



< v < 1 
^ > 1 

< i/ < 1 
v > 1 



(23) 



(24) 



(We note that B u {2, 1) = For Dotri scenarios, we had V£ t s = = giving 

us P^p = jq — which is the smallest probability we will report in this entire paper. 

So, in our first set of simple (m = 5) scenarios, we observe a decrease in the probabilities 
of separability from the real to the complex to the quaternionic well as a progression 

from ^ to v to v 2 in the functional forms occurring in the corresponding HS separability 
probability functions. 

The exponents of v in this progression, that is |,1,2 bear an evident relation to the 



Dyson indices 34 



2, 4, corresponding to the Gaussian orthogonal, unitary and sym- 



plectic ensembles 35]. Certainly, this observation bears further investigation, in particular 



since the foundational work of Zyczkowski and Sommers [17|) in computing the HS (sepa- 
rable plus nonseparable) volumes itself relies strongly on random matrix theory. ("[T]hese 
explicit results may be applied for estimation of the volume of the set of entangled [emphasis 
added] states" 17|, p. 10125].) However, this theory is framed in terms of the eigenvalues 
and eigenvectors of random matrices — which do not appear explicitly in the Bloore pa- 
rameterization — so, it is not altogether transparent in what manner one might proceed. 
(But for the m = 5 highly sparse density matrices for this set of scenarios, one can explicitly 
transform between the eigenvalues and the Bloore parameters.) 
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B. Four nullified pairs of off-diagonal entries — 15 scenarios 



1. 5- dimensional real case — PS.„ = I; t4s 



Here, there are fifteen possible scenarios, all with V£ t s = ^q. Six of them are trivial 
(separability probabilities of 1), in which c scenario is either n (scenarios [(1,2), (1,3)], [(1,2), 
(2,4)], [(1,3), (3,4)] and [(2,4), (3,4)]) or 4 (scenarios [(1,2), (3,4)] and [(1,3), (2,4)]). Eight of 
the nine non-trivial scenarios all have — similarly to the 4-dimensional analyses (sec. Ill A lj) 
— - separability functions <S(z/) either of the form, 



S scenarioiy^) — \ ; (^5) 




(for scenarios [(1,2), (2,3)], [(1,3), (2,3)], [(2,3), (2,4)] and [(2,3), (3,4)]) or, dually, 

{TV < V < 1 
(26) 

(for scenarios [(1,2), (1,4)], [(1,3), (1,4)], [(1,4), (2,4)] and [(1,4), (3,4)]). The corresponding 
HS separability probabilities, for all eight of these non-trivial scenarios, are equal to | = 
0.625. This result was, in all the eight cases, computed by taking the the ratio of V^f = ^ 



lu v tot — 480- 

In the remaining (ninth) non-trivially entangled case — based on the non-nullified dyad 
[(1,4), (2, 3)] — we have, taking the ratio of Vj^f = ^ to V t ^ t s = a quite different Hilbert- 
Schmidt separability probability of ^ ~ 0.54038. This isolated scenario (with c scenario = 4) 
can also be distinguished from the other eight partially entangled scenarios, in that it is the 
only one for which entanglement occurs for both v < 1 and v > 1. We have 




■%l,4),(2,3)](f) = < • (27) 



By way of illustration, in this specific case, we have the scenario-specific marginal jacobian 
function, 

q ( , 07(-3z/ 2 + (^ + 4) + l)log(/y)+3) 

J[(1,4),(2,3)]W = 30(1/ - l) 5 ' " ' 



14 



2. 6- dimensional mixed (real and complex) case — P t 



HS _ 105tt . 135tt . 3 
sep 512 ' 1024 ' 8 



Here, we again nullify all but two of the off-diagonal entries (m = 4) of p, but allow the 
first of the two non-nullified entries to be complex in nature. Making (apparently necessary) 
use of the circular /trigonometric transformation pu = r 2 sin# 2 , P22 = f 2 cos# 2 , we were able 
to obtain an interesting variety of exact results. One of these takes the form, 



S [(&)M)M = 5 [(i;3),(l,4)]M = {{f ,0 < V < 1} {$-,V > 1} . 



(29) 



Now, we have V t g s = ^ and V s % s = Jgj, so P s % s = f£ « 0.644272. The two dual 
scenarios — having the same three results — are [(1, 2), (2, 3)] and [(1, 3), (2, 3)]. 
Additionally, we have an isolated scenario, 

5 [(i;4),(2,3)](^) = {{2^0<^<1} {f,u>l}, (30) 

for which, V* t s = ^ and V s f/ = Jf^, so P£/ = ±§f « 0.414175. (Note the presence of 
&o£/j and z/ in (!30|) — apparently related to the mixed [real and complex] nature of this 
scenario (cf. fl33l ).) 
Further. 

5 [(i;4) >( 2,4)](^) = %4),(3,4)](^) = {{f ,0 < 1/ < 1} {£,1/ > l} , (31) 

the dual scenarios being [(2, 3), (2, 4)] and [(2, 3), (3, 4)]. For all four of these scenarios, 

V HS — 71-2 anrl V HS — 7f2 cn P HS — 5 _ n ^7^, 
v tot — 1440 dI1U V sep — 3840' bU sep — 8 ~~ 



5. 7- dimensional complex case — -P S ep = g 

Here, in an m = 4 setting, we nullify four of the six off-diagonal pairs of the 4x4 density 
matrix, allowing the remaining two pairs both to be complex. We have (again observing a 
shift from ^/v in the real case to v in the complex case) 

J £ < ^ < 1 

< ^[(1>),(1~4)]( I/ ) = ^[(MMM)]^) = < ^[(1~4),(2>)]( Z/ ) = ^[(1~4),(3~4)] I 17 ) = ^ „ • (3 2 ) 

U- 1/ > 1 
I 2i/ 

Since V t ^ t s = and V^f = j^q, we have Pj™ = | = 0.4. We have the same three out- 
comes for the four dual scenarios [(1~2), (2~3)], [(1~3), (2~3)], [(2~3), (2~4)] and [(2~3), (3~4)], 
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as well as — rather remarkably — for the (again isolated [cf. (j3"01) ]) scenario [(1,4), (2,3)], 
having the (somewhat different) separability function (manifesting entanglement for both 
v < 1 and v > 1), 



TT 2 U < V < 1 

s \a~A\(za)]( u ) = { , ■ ( 33 ) 



(However, c scenario = tt 2 for this isolated scenario, while it equals ^ for the other eight.) 
The remaining six (fully separable) scenarios (of the fifteen possible) simply have = 1. 

4- 8- dimensional mixed (real and quaternionic) case 
We report here that 

C [(I^),(l,4)] = "Tg"' C [(1,2),(M)] = 32 ' ( 34 ) 

where as before the wide tilde notation denote the quaternionic off-diagonal entry. 



C. Three nullified pairs of off-diagonal entries — 20 scenarios 

1. 6- dimensional real case — P^ = 2 — |^|| ; rj ^ 

Here (m = 3), there are twenty possible scenarios — nullifying triads of off-diagonal pairs 
in p. Of these twenty, there are four totally separable scenarios — corresponding to the 
non-nullified triads [(1,2), (1,3), (2,4)], [(1,2), (1,3), (3,4)], [(1,2), (2,4), (3,4)] and [(1,3), 
(2,4), (3,4)] - - with c scenario = ^ and Vgf = V™ = To proceed further in this 

6-dimensional case — in which we began to encounter some computational difficulties - 
we sought, again, to enforce the four nonnegativity conditions (<^, (J1J, (j3j)), but only 
after setting v to specific values, rather than allowing v to vary. We chose the nine values 
v = |,|,|,|, 1, 2, 3, 4 and 5. Two of the scenarios (with the triads [(1,2), (2,3), (3,4)] 
and [(1,3), (2, 3), (2, 4)]) could, then, be seen to fit unequivocally into our earlier observed 
predominant pattern, having the piecewise separability function, 

fzE^E < v< 1 

5[(l,2),(2,3),(3,4)](l / ) = S[(l,3),(2,3),(2,4)](^) = < ■ (35) 

It 
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FIG. 1: Plot of iSrn ) 2),(i,4),(2,3)] ( v ) an d the close fitting function ([37]) . The former function slightly 

2 

dominates the latter except at v = 1, where they both equal ^- « 4.9348. 

We, then, computed for these two scenarios that V£ t s = gfgg « 0.00538303 and (again mak- 
ing use of the transformation p n = r 2 sin 6 12 , p 2 2 = r 2 cos # 2 ) that V^f = 2 ^gfgQ — 736432 
0.00358207. This gives us P S H J = 2 - f||f » 0.665437. For two dual dyads, we have the 
same volumes and separability probability and, now, the piecewise separability function, 

f < 1/ < 1 

«%,2),(l,4),(3,4)](y) = «5[(1,3),(1,4),(2,4)] (^) = ^ ■ (36) 

Additionally, in Fig. [TJ we are able to plot «S[(i j2 ) ) (i J 4) ) (2,3)](i / ) along with its close fit to 

^ < v< 1 
S m {y) = { . (37) 

(The analogous plots for the scenarios [(1,3), (1,4), (2,3)] and [(1,2), (1,3), (2,3)] appear 
to be precisely the same in character as Fig. [TJ) If we use the close fit (|37|) as a proxy 
for «S[(i )2 ) ) (i,4) j (2,3)](i / ), we obtain an approximate HS separability probability of ^ = J^ 2 - = 

5760 

0.5625. 

We have not, to this point, been able to explicitly and succinctly characterize the functions 
<S 'scenario^) for non-trivial fully real m = 3 scenarios other than the dual pair ( fl35l) . (136])). 

In all the separably non-trivial scenarios so far presented and discussed, we have had the 
relationship S scenario {l) = c scenario . However, in our present m = 3 setting (three pairs of 
nullified off-diagonal entries), we have situations in which S scenar i {l) < c scenar i . The values 
of c scen ario in the sixteen non-trivial fully real m = 3 scenarios are either ^ ~ 4.9348 (twelve 
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occurrences) or f « 4.18879 (four occurrences — [(1,2), (1,3), (1,4)], [(1,2), (2,3), (2,4)], 
[(1,3), (2,3), (3,4)] and [(1,4), (2,4), (3,4)]). In all four of the latter occurrences, though, 
we have the inequality, 

1 47T 
SscenarM = ^ ( 12 + 16?r + ^) « 3 ' 8281 < "J « 4.18879, (38) 

2 

as well as a parallel inequality for four of the twelve former (^-) cases. The implication 
of these inequalities for those eight scenarios is that at v — 1 (the value associated with 
the fully mixed [separable] classical state), that is, when pnp^ = P22P33, there do exist 
non-separable states. 



2. 7- dimensional mixed (one complex and two real) case — Pj^jj = ig 

Here, in an m = 3 setting, we take the first entry of the non-nullified triad 
to be complex and the other two real. Of the twenty possible scenarios, four — - 
[(1~2), (1, 3), (1, 4)], [(1~2), (2, 3), (2, 4)], [(1~3), (2, 3), (3, 4)] and [(1~4), (2,4), (3,4)] - - had 

2 

Cscenario — \ ~ 4.9348 and these four all had the same (lesser) value of 

^6 7T 2 

S scenario (l) = — + — « 4.54148. (39) 

There were seven scenarios with c scenar io = -f 2 - ~ 5.58505. Three 
of them [(1,2), (1,3), (2,4)], [(1,2), (1,4), (2,3)] and [(1,3), (1,4), (2,3)] 

had S scenaT i {l) = ^ (manifesting equality), while four 
[(1,2), (1,3), (2,3)], [(^2), (1,4), (2,4)], [(1,3), (1,4), (3,4)] and [(2~3), (2,4), (3,4)] 
had the result (|39|) (manifesting inequality). 

The remaining nine of the twenty scenarios all had c scenaT io = S scenaT i {l) = ^j- ~ 6.57974. 
For one of them, we obtained 

g [(i;2),(2,3),(3,4)](^)=| <I/<l' (40) 

with associated values of V^ 5 = t^o' = 322560 an d = i§ ~ 0.6875. A dual 

scenario to this one that we were able to find was [(1, 2), (1, 4), (3, 4)]. The separability func- 
tions — and, hence, separability probabilities — for the other eighteen scenarios, however, 
are unknown to us at present. 
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3. 8- dimensional mixed (two complex and one real) case 
Our sole result in this category is 

c [(i;2),(i; 3 ),(i,4)] = - jQ-- ( 41 ) 

4- 9-dimensional complex case 

Now, we have three off-diagonal complex entries, requiring six parameters for their spec- 
ification. According to the remarks of Strzebonski (sec. II Cp . this is about the limit in the 
number of free off-diagonal parameters for which we might hopefully be able to determine 
associated separability functions. 

As initial findings, we obtained 

i ^[(i;4),(2;3),(2;4)]( 1 ) = c [(i;4),(2;3),(2;4)] = ~^-> (42) 

and also for scenarios [(1~2), (1~4), (3~4)], [(1~3), (1~4), (2~4)] and [(1~4), (2~3), (3~4)], while 

TT 3 

c [(i^),(i^),(M)] = c [(i;3),(2; 3 ),(3;4)] = c [(i;4),(2; 3 ),(2;4)] = y • ( 43 ) 

D. Two or fewer nullified pairs of off-diagonal entries 

1. 7- dimensional real case 

The [(1,2), (1,3), (2,4), (3,4)] scenario is the only fully separable one of the fifteen possible 
(m = 2). For all the other fourteen non-trivial scenarios, there are non-separable states both 
for v < 1 and v > 1. For all fifteen scenarios, we have c scenario = 2|i pa 6.57974. Otherwise, 
we have not so far been able to extend the analyses above to this m = 2 fully real case (and 
a fortiori the m = 1 fully real case), even to determine specific values of S scenario {l). 

2. 8- dimensional real case 

Here we have c scenaT io = ^f- ~ 8.77298 for all the six possible (separably non-trivial) 
scenarios (m = 1). Let us note that this is, in terms of preceding values of these constants 
(for the successively lower-dimensional fully real scenarios), = f(^f-), while ^ = 
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Also, = |(^|-), the further relevance of which will be apparent in relation to our 



discussion of the full 9-dimensional real scenario (sec. IVIIIj) 



III. QUBIT-QUTRIT ANALYSES 

The cancellation property, we exploited above, of the Bloore parameterization — by which 
the determinant and principal minors of density matrices can be factored into products of 
(nonnegative) diagonal entries and terms just involving off-diagonal parameters (zij) - 
clearly extends to n x n density matrices. It initially appeared to us that the advantage of 
the parameterization in studying the two-qubit HS separability probability question would 
diminish if one were to examine the two-qubit separability problem for other (possibly 
monotone) metrics than the HS one, or even the qubit-qutrit HS separability probability 
question. But upon some further analysis, we have found that the nonnegativity condition 
for the determinant of the partial transpose of a real 6x6 (qubit-qutrit) density matrix (cf. 
(J2])) can be expressed in terms of the corresponding z^-'s and two ratio variables (thus, not 
requiring the five independent diagonal variables individually), 

v x = W™, v 2 = ^» (44) 
P22P44 P33P55 

rather than simply one (y) as in the 4x4 case. (We compute the qubit-qutrit partial 
transpose by transposing in place the four 3x3 blocks of p, rather than — as we might 
alternatively have done — the nine 2x2 blocks.) 

A. Fourteen nullified pairs of off-diagonal entries — 15 scenarios 

1. 6- dimensional real case — = || 

To begin our examination of the qubit-qutrit case, we study the (m = 14) scenarios, in 
which only a single pair of real entries is left intact and all other off-diagonal pairs of the 
6x6 density matrix are nullified. (We not only require that the determinant of the partial 



transpose of p be nonnegative for separabi 



ityto hold — as suffices in the qubit-qubit case, 



given that p itself is a density matrix 36j, |37| - but also, per the Sylvester criterion, a 



nested series of principal leading minors of p.) We have six separably non-trivial scenarios. 

7T 

1440 - 
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(For all of them, Vg s = fig.) 



Firstly, we have the separability function, 



The dual scenario to this is [(2,4)]. Further, 




v x > 1 



V\ < 1 



(45) 




^1^2 < 1 



V\V 2 > 1 



(46) 



with the dual scenario here being [(3,4)]. Finally, 




2 



2 



v 2 > 1 



v 2 < 1 



(47) 



having the dual [(3,5)]. 

The remaining nine possible scenarios — the same as their complex counterparts in the 
immediate next analysis — are all fully separable in character. 

We have found that V s ^ p s = for the six non-trivially separable scenarios here, so 
P^p — f| ~ 0.589049, as in the qubit-qubit analogous case (sec. Ill A II) . 

2. 7 -dimensional complex case — P^p = | 

Now, we allow the single non-nullified pair of off-diagonal entries to be complex in nature 
(the two paired entries, of course, being complex conjugates of one another). (V t ^ t s = -g^j 
for this series of fifteenb scenarios.) Then, we have (its dual being [(2,4)]) 




(48) 



Further, we have (with the dual [(3,4)]) 




(49) 



and (its dual being [(3,5)]) 




(50) 
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For all six of these scenarios, V sep = j^o, so P^s — §• 

B. Thirteen nullified pairs of off-diagonal entries — 105 scenarios 

1. 7- dimensional real case — Pj^j = |; ^; ||; 

Continuing along similar lines (m = 13), we have 105 combinatorially distinct possible 
scenarios. Among the separably non-trivial scenarios, we have 



TV V\ < 1 

77= 



5 [(f2)il.5)l(^)- 5 [ff4Ul,5)l( i/ l)-^ ' ( 51 ) 



(dnals being [(1,2),(2,4)] and [(1,4),(2,4)]). We computed V* 8 = ^,V S % S = so 
Psep = I = 0-625 for these scenarios. 
Also, 



5 r 8 u,3Ui,6ii V2 ) - 5 f(S),(i,6)i (^i' "2) - < • (52) 



7T Z/xZ/ 2 < 1 

Z/lZ/ 2 > 1 



We, then, have V* s = K?/ = gfe, so Pf* = £ = 0.3125. 
Additionally, 



4 i/ 2 < 1 

4, 



^1.4^1-2.6)1 (^2) - < • (53) 



'[(1,4),(2,6)] 

For this scenario, we have l^ 5 = , V™ = ^j, so Pf* = § » 0.294524. 
Further still, 

[ ^1 < 1 

For this scenario, we have V« t s = V^/ = 3^, so = $ » 0.54038. 
Further, 



- 7T 1^2 < 1 

5 [a 6 2),(2,6)](^) = ^ . . • (55) 

, 2 ( cos (v T ^) + ^ 1 ) " 2>1 



The separability function for [(1,3), (1,5)] is obtained from this one by replacing v 2 by v\. 
Also, 

5 KiW)](^.^)= (56) 
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4 



^1^2 < 1 
2/^2 > 1 



The separability function for [(1,2), (3, 5)] can be obtained from this one by setting v x = 1. 



2. 8- dimensional mixed (real and complex) case - 



HS _ 105tt 
sep 512 



Further, we have (with V t 



HS _ 7T 2 



tot 80640 



for all scenarios) 



S 6X - 6 (^) = s 6x ~ 6 {Vl 

V U [(1,4), (2,4)] V X ' 



[(1.2),(2,4)]' 



f 

^ 0<^<1 



Since V^f = ggf^g, we have P s ^f = « 0.644272 for both these scenarios. 



Further, 



-6x6 



^ < v x < 1 



and 



-6x6 



[(1,3), (3,4)] V lj 2J [(1,4), (3,4)] V 11 11 



3 



V\^2 > 1 



Additionally, 



i- 

3 



\iry/vxv 2 < v x v 2 < 1 



VlV 2 > 1 



\-Ky/ViV 2 (3 - ^1^2) < Z/1^2 < 1 



and 



3 



V\V 2 = 1 

< V X V 2 < 1 • 



We have also obtained the separability function (Fig. [2]) 

v x > 1 



?6x6 



[(2>),(2,4)]V *' 



4-7T 

3 



§7T (3 - V X ) y/Vx <Vx <1 



Of the 105 possible scenarios, sixty had <S^ C g^ aHo (l, 1) 



(57) 



(58) 



(59) 



(60) 



(61) 



(62) 



Cscenario = ^, thirty-three had 



6x6 



scenario 



4tt 
3 



: C 

2tt 



cenamo 



2n, and twelve (for example, [(3,4), (5,6)]) had S 



6x6 /-1 i> 
scenarioK 1 i 
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FIG. 2: Plot of the separability function S" 6 *-? , (u^ 

J [(2,3), (2,4)] V U 



3. 9- dimensional complex case — P^ijj = |; | 
We have obtained the results 



- > 1 

Since V^ t s = and V™ = we have here P S H J = § = 0.4. We have the same 

three outcomes also based on the separability function, 



, .r ^1^2 > 1 



Further, 



and 



-t{ 1 v 1 v 2 < Z/iZ/ 2 < 1 



2 



< v x < 1 



- 2 -2^ 



For both of these last two scenarios, we have V^ 5 = 36 ^ 88Q and V^f = 36 2 880 , leading to 
Psep = I ~ 0.33333. Also, we have these same three outcomes based on the separability 
function, 



27rViZ/ 2 (2 - ^l^ 2 ) < V X V 2 < 1 
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Of the 105 possible scenarios — in complete parallel to those in the immediately pre- 
ceding section — sixty had S^ ario (l, 1) = c scenario = ^, thirty-three had S^ ario (l, 1) = 
Cscenario = ^ ! , and twelve (for example, [(3, 4), (5, 6)]) had S^ e l ario (l, 1) = ^ < c scenario = 

7T 2 . 

Our results in this (9-dimensional) section and the (8- dimensional) one immediately pre- 
ceding it are still incomplete with respect to various scenario-specific separability functions 
and, thus, the associated HS separability properties. 



C. Twelve nullified pairs of off-diagonal entries — 455 scenarios 



1. 8- dimensional real case 



Now, we allow three of the off-diagonal pairs of entries to be non-zero, but also require 
them to be simply real. We found the separability function 



«5[(l,2),(l,3),(3,4)](^l, v 2) = <%1,2),(1,4),(3,4)](^1, V ?) 



(68) 



3 



:7T 



^ |tt (3 sec- 1 (y/w) + A^U~ 2 + ^El _ i K /T7^ | 
Also, we have 



<%,3),(1,4),(2,4)] (^l) = < 



-Itt 
3 

3 



vr (4i/ 1 3/2 + (3 sin- 1 [J^~] -Mfcpt) v x + v / ^l) 



3i/i 



J- = !/ 2 A z/i > 

v x > A ^ > f 2 A f 2 > • 

z/i > A i < i/ 2 

z/x = 1 
< vx < 1 ■ 
v x > 1 



(69) 



IV. QUTRIT-QUTRIT ANALYSES 



In the qubit-qubit (4x4 density matrix) case, we were able to express the condition (JHJ) 
that the determinant of the partial transpose of p be nonnegative in terms of one supplemen- 
tary variable (v), given by ([6]), rather than three independent diagonal entries. Similarly, in 
the qubit-qutrit (6x6 density matrix) case, we could employ two supplementary variables 
(vi,v 2 ), given by (1441) . rather than five independent diagonal entries. 
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For the qutrit-qutrit (9x9 density matrix) case, rather than eight independent diagonal 

entries, we found that one can employ the four supplementary variables, 

P11P55 P22P66 P44P88 P55P99 ,„ n \ 
v\ = ; v 2 = ; V3 = ; ^4 = • (70) 

P22P44 P33P55 P55P77 P66P88 

A. Thirty-five nullified pairs of off-diagonal entries — 36 scenarios 

1. 10- dimensional complex case — Pppx = \'i \ 

Here, we nullify all but one of the thirty-six pairs of off-diagonal entries of the 9x9 
density matrix p. We allow this solitary pair to be composed of complex conjugates. Since 
the Peres-Horodecki positive partial transposition (PPT) criterion is not sufficient to ensure 
separability, we accordingly modify our notation. 

Our first result is 



, TT V\ < 1 

s S„<">> = " < 71 > 

1 ^ v x > 1 



(a dual scenario being [(2,4)]). We have V t * s = g^fgoo, V^ 5 T = los £ mo , so P^ T = §. 
The same three outcomes are obtained based on the PPT function 

{7T < 1 

. <72) 
— v x v 2 > 1 

On the other hand, we have V t ^ t s = 362 g 80 o 1 Vppt = 21772800 ' an< ^ ^ppt = | based on the 
PPT function 



,71 Va > 1 

S?*_ 9 (uA ={ " . (73) 

[(6,8)] V I V ' 

7TZ/ 4 < 1/4 < 1 



Of the thirty-six combinatorially possible scenarios, thirteen had Ppp T = |, while four 
had Ppp T = |, and the remaining nineteen were fully separable in nature. 



B. Thirty- four nullified pairs of off-diagonal entries — 630 scenarios 

1. 12- dimensional complex case — Pppx = 3> M 

Since the number of combinatorially possible scenarios was so large, we randomly gener- 
ated scenarios to examine. 
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Firstly, we found 



tH u 2 > 1 



n z u 2 < v 2 < 1 

For this scenario, we had V t % s = 479 ^ W00 , V^ T = 1437 g 4800 , giving us P#/ T = §. 
Also, we found 

z/ 2 z/ 4 < 1 

^ : , • (75) 

For this scenario, we had V t % s = 479 ^ 1600 , Vg? T = ^filooo' S ivin g us P ppt = ^ ~ 0.23333. 



V. QUBIT-QUBIT-QUBIT ANALYSES, I 

For initial relative simplicity, let us regard an 8 x 8 density matrix p as a bipartite system, 
a composite of a four-level system and a two-level system. Then, we can compute the partial 
transposition of p, transposing in place its four 4x4 blocks. The nonnegativity of this partial 
transpose can be expressed using just three ratio variables, 

P11P66 P22P77 P33P88 , nrx 

v\ = ; vi = ; ^3 = , (76) 

P22P55 P33P66 P44P77 

rather than seven independent diagonal entries. 

A. Twenty-seven nullified pairs of off-diagonal entries — 28 scenarios 

1. 9- dimensional complex case — PppT = \ 
We have the PPT function 



(Scenario [(2,5)] was dual to this one.) For this scenario, V t ^ t s = 36 2 880 , Vpp T = 108 g 640 > 
yielding Ppp T = |. There were twelve scenarios, in toto, with precisely these three outcomes. 
The other sixteen were all fully separable in nature. 
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B. Twenty-six nullified pairs of off-diagonal entries 



— 378 scenarios 



1. 11- dimensional complex case — Pppj ~ ' " ' 



ppt — 3 ' 9 



Again, because 
inspection. 

Firstly, we obtained 



of the large number of possible scenarios, we chose them randomly for 



5 [(3; 5 ),(6; 8 )]^^ 2 ) - i 2 nA i A n 

i\ l v x v 2 v x > A — > v 2 A v 2 > 



(78) 



V - 1 

(Of course, the symbols "A" and "V", used by Mathematica in its output, denote the logical 
connectives "and" (conjunction) and "or" (intersection) of propositions.) For this scenario, 

. __7rn _2 __7-r^ _2 . . _jrn 1 



we had V^f = 
Also, 



___ v HS = — — 

39916800' v PPT 119750400 



giving us Ppp T = |. 



7T 2 


^i > 1 A u 3 > 1 




0<//i<lA^ 3 >l 


7T 2 Z/ 3 


f i > 1 A < y 3 < 1 


7T 2 Z/il/ 3 


0<i/i<lA0<y 3 <l 



(79) 



For this scenario, we had V£ t s = 3g9 f emo , V^ T = 3 ^ 1200 , giving us P^ T = §. 
We also found the PPT function 



[(1,3),(4,7)] V 1J ^ " 



u 3 = 1 

< f 3 < 1 



(80) 



7 ( C ° S 1 (V 1- A) ~ CSC 1 (v^O) ^3 



VI. QUBIT-QUBIT-QUBIT ANALYSES. II 



Here we regard the 8x8 density matrix as a tripartite composite of three two-level 
systems, and compute the partial transpose by transposing in place the eight 2x2 blocks 
of p. (For symmetric states of three qubits, positivity of the partial transpose is sufficient 
to ensure separability jss]].) Again the nonnegativity of the determinant could be expressed 
using three (different) ratio variables, 



PnPu, 
P22P33 



PUP5B 
P33P66 

28 



^3 



P55P88 

PmPn 



(81) 



A. Twenty-seven nullified pairs of off-diagonal entries — 28 scenarios 



1. 9- dimensional complex case — Pppt = \ 



There were, again, twelve of twenty-eight scenarios with non-trivial separability proper- 
ties, all with V* t s = Vfigr = logfeio' yielding P^ T = §. One of these was 



71 V\ < 1 



^ V X > 1 



(82) 



2. 11- dimensional complex case — PppT = m">k 



We obtained the PPT function 



7T 



41^1^3 



^ = U3 A z/ x > A v 2 > 

zv 2 > A ((z/i = A u 3 > 0) V (z/ 3 = A v x > 0)) 
z/ x > A v 2 > A £i < z/ 3 



For this we had V^ s = 39^00, ^/pt = 



^2 _ n^s ^ > A z/ 2 > A ^ > z/ 3 A z/ 3 > 



17tt 2 



239500800 



Additionally, 



<S 8x8 - (i/x) 

[(M),(W U 



7T 



y\ < 1 

V 1 >1. 



Here, we had V^ 5 



^pt = 1197^00 ' giving us P^/ T 



tot 39916800 ' 

Another PPT function we were able to find was 



(83) 



, giving us P^ T = g « 0.283333. 



(84) 



7T 

2 



i>2 > A (1/3 = z/ 2 V (z/ 2 > ^3 A v 2 < 2z/ 3 ) V {y 2 > 2z/ 3 A z/ 3 > 



^gf^ "2>0Av 2 <v 3 



(85) 
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VII. APPROXIMATE APPROACHES TO 9-DIMENSIONAL REAL QUBIT- 
QUBIT SCENARIO 



As we have earlier emphasized (sec. II Cp . it appears that the simultaneous computational 
enforcement of the four conditions ((|2J), Q, ©) that would yield us the 9-dimensional 
volume of the separable real two-qubit states appears presently highly intractable. But if 
we replace (|5]) by less strong conditions on the nonnegativity of the partial transpose {p T ), 
we can achieve some form of approximation to the desired results. So, replacing (j3J) by the 
requirement (derived from a 2 x 2 principal minor of p T ) that 

1 - vz\ A > 0, (86) 

we obtain the approximate separability function (Fig. |3J) 



5127T 2 
27 



< v < 1 



SrealW) = \ , , ,x ■ (87) 

27^ ^ ">1 



(In the analyses in this section, we utilize the integration limits on the z^s [15l . eqs. (3)-(5)] 
yielded by the cylinrical decomposition algorithm [CAD], to reduce the dimensionalities of 
our constrained integrations.) This yields an upper bound on the separability probability 
of the real 9-dimensional qubit-qubit states of \ + 1 ^ 2 ~ 0.88427. We obtain the same 
probability if we employ instead of (186]) the requirement 



v - 4 > 0, (88) 
which yields the dual function to (187]) . namely, 



512tH 



, - V > 1 

Sreal{v) = { • (89) 

-^7r 2 (z/-3)v^ < i/ < 1 

(The left-hand sides of f)86p and f)88p are the only two of the six 2x2 principal minors 
of p T that are non-trivially distinct — apart from cancellable nonnegative factors — from 
the corresponding minors of p itself.) If we form a "quasi-separability" function by piecing 
together the non-constant segments of (1561) and (155]) . we can infer — using a simple sym- 
metry, duality argument — an improved lower bound on the HS separability probability of 
w 0.76854. 
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Sreal (v) 




FIG. 3: Approximation (|87p to the separability function for the 9-dimensional real two-qubit states 



A. Analytically-derived beta functions 



We have not, to this point of time, been able to jointly enforce the constraints (1561) 
and (1551) . (Doing so, would amount to solving a constrained /owr-dimensional integration 
problem.) However, we had some interesting success in jointly enforcing (1561) along with the 
constraint 

- zf 2 + 2^/uz 1 ^z li z 12 - zl 3 - vz 2 u + 1 > 0, (90) 

corresponding to the first of the four principal 3x3 minors of p T . (This amounts to solving 
a constrained t/iree-dimensional integration problem.) In attempting to jointly enforce these 
two constraints, we obtained (without, however, being directly able to explicitly derive the 
separability function over /i e [0,oo]) an upper bound on V^j? (and a consequent upper 



bound on of 0.77213) of 

llvr 4 



V| , ronn + / Jreal(l^) Sreal ~ 0.00124359. 

zllooUU ./q 



(91) 



Here we have the approximate separability function (monotonically-c/ecreasmgfor \i e [0, 1]), 



Sreal 



7T 



(135 M2 (I |) /i 2 + (2^ 2 - 13) (|, f) - 2B,i (|, I)) 



3/x 3 



(92) 



For (apparently necessary) convenience in computation, so that Mathematica only has to 



deal with integra 
original study 



iJreal (aO 



powers rather than half-integer ones, we use fi = ^fv (as we did in our 
27]]) and the corresponding jacobian function 0, eq. (10)] (cf. (115]) ) is 
/i 4 (12 ((n 2 + 2) (// + 1 V + 8) fi 2 + 1) log(/i) - 5 (5fi s + 32/i 6 - 32^t 2 - 5)) 



1890 (p 2 - V 



(93) 
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What is most interesting here, of course, is that we have beta functions (now derived ana- 
lytically) explicitly appearing in the (approximate) separability function formula (I921) . while 



we only previously had numerical evidence for their relevance [15 ]. 

Though we sought to enforce additional constraints based on the principal minors of 
p T , in order to obtain tighter upper bounds on V^f and -Pjfjf, this did not seem to be 
computationally possible. (One might also pursue similar "approximate" strategies in the 
lower-dimensional instances, studied earlier in this paper, still not apparently amenable to 
exact solutions.) 

Also, as was mentioned in the Introduction, we have not yet been able to exploit to derive 
HS separability functions and probabilities, the simplified integration limits fill I) based on 
the spheroidal parameterization ffTUj) suggested by the geometric discussion of Bloore 
sees. 6, 7]. 

VIII. CONCLUDING REMARKS 

Though we have regarded the derivation of the scenario-specific separability (and PPT) 
functions as our primary theoretical and computational challenge, the derivation of the HS 
separability (and PPT) probabilities from such functions — that is, the performance of an 
integration over the (n-l)-dimensional simplex spanned by the diagonal entries of the n x n 
density matrix p — can also be quite difficult. (As a general, somewhat informal observation, 
real off-diagonal entries in our scenarios tend to yield square roots of the diagonal entries in 
the integrands, which prove more challenging to integrate than complex off-diagonal entries.) 
The qubit-qubit results above motivated us to reexamine previously obtained results (cf. 



271 . eqs. (12), (13)]) and we would like to make the following observations pertaining to the 
full 9-dimensional real and 15-dimensional HS separability probability issue. We have the 
exact results in these two cases that 

rco rl 2 

J Jrealiy) = 2 J J real (u) = um8Q « 8.60561 • 1(T 6 (94) 

and 

JcompUv) = 2 J JcompUv) = « 9 - 91072 ' . ( 95 ) 

Now, to obtain the corresponding total (separable plus nonseparable) HS volumes 171 ]. that 
is, » 0.0016106 and 8513 f 0500 « 1.12925 ■ 10~ 6 , one must multiply (El) and (J95D by the 
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factors of C r , 



real 



512tH 



2V 



-i 7 ~~tt" ~ 187.157 and C comp i ex ^7 — —rr 
To most effectively compare these previously-reported results with those derived above 
in this paper, one needs to multiply C rea i and S rea i(u), by 2~ 4 = i and in the complex case 



2 7T 



1139.42, respectively. 



by 2" 



io. Doing so, for example, would adjust C rea i to equal c, 



128 



real 



32tH 
27 



11.6973, 



which we note, in line with our previous series of calculations [sec. Ill D 2j is equal to |(^f-)- 
(Andai 39J also computed the same volumes — up to a normalization factor — as Sommers 
and Zyczkowski [l^.) Now, our estimates from 15] are that S rea i(l) = 114.62351 < C rea i 
and S comp i ex {l) = 387.50809 < C comp i ex . These results would appear — as remarked above 
- to be a reflection of the phenomena that there are non-separable states for both the 9- 
and 15-dimensional scenarios at v — 1 (the locus of the fully mixed, classical state). 
If we take trial function, 



S real (u) = 114.62351 



V" < v < 1 



\ 



4= v > 1 



(96) 



/ 



and insert it into our formula f|T3|) . as at least the lower-dimensional real results in Fig. [T] 
and (1271) might suggest trying, we obtain the outcome, V s ^ real = 0.000613694 and (using 
the known V^ reoX 



171 ]) an associated HS separability probability of 0.381034. This — as we 



expected — is certainly not in accord with the numerical results of 1_5|, sec. V.A.2]. The HS 
real separability probability was estimated there to be P^ ep / rea i = 0.4538838 and S rea i(v) (cf. 
(fl2"]l ) was estimated to be closely proportional to the incomplete beta function B v (^, a/3) 
and not the quite different function y/u oc B v (^, 1), as our small exercise here assumes. 

So, we may say, in partial summary that we have been able to obtain certain exact 
two-qubit HS separability probabilities in dimensions seven or less, making use of the ad- 
vantageous Bloore parameterization Q, but not yet in dimensions greater than seven. This, 
however, is considerably greater than simply the three dimensions (parameters) we were able 
to achieve [3] in a somewhat comparable study based on the generalized Bloch representa- 
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4l|. In 







extending an approach of Jakobczyk and Siennicki 



tion parameterization 

4l| - we primarily studied two- dimensional sections of a set of generalized Bloch vectors 
corresponding tonxn density matrices, for n = 4, 6, 8, 9 and 10. For n > 4, by far the most 
frequently recorded HS separability [or positive partial transpose (PPT) for n > 6] proba- 
bility was j ~ 0.785398. A very wide range of exact HS separability and PPT probabilities 
were tabulated. 
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Immediately below is just one of many matrix tables (this one being numbered (5) in 
Q]) presented in 7] (which due to its large size has been left simply as a preprint, rather 
than submitted directly to a journal). This table gives the HS separability probabilities 
for the qubit-qutrit case. In the first column are given the identifying numbers of a pair 
of generalized Gell-mann matrices (generators of SU(6)). In the second column of (1971) are 
shown the number of distinct unordered pairs of SU (6) generators which share the same 
total (separable and nonseparable) HS volume, as well as the same separable HS volume, 
and consequently, identical HS separability probabilities. The third column gives us these 
HS total volumes, the fourth column, the HS separability probabilities and the last (fifth) 
column, numerical approximations to the exact probabilities (which, of course, we see - 
being probabilities — do not exceed the value 1). (Due to space/page width constraints, we 
were unable to generally present in these data arrays the HS separable volumes too, though 
they can, of course, be deduced from the total volume and the separability probability.) 



I {1,13} 48 

{3,11} 4 

{3,13} 4 

{3,25} 4 

{8,13} 4 

{8,25} 4 

{11,15} 4 

{11,24} 2 

{13,24} 2 

{13,35} 4 

{15,16} 4 

{16,24} 2 

{20,24} 2 

{24,25} 2 

{24,27} 2 

I {25,35} 4 



i 

9 

27 

4 
9 

27 

2 

3 

V2 
3 

AV2- 



7T 



27 
2Vf 
72 



25 



72 



5 

32^ 
81 

25 
144 

25 
144 

25 
27^ 

25 
27^ 

V3tt 
5 



7T 

4 

1 

V2 
5 
6 

5 

AV2 
1 

V3 
l~2 

3 

l , 3yjs 

3 ' 4tt 

f + §sin 



_8_ 

75 



Hi) 

(-2 + 5^) 



-i (j» 



TV 



7T 



Y2 (5 + 3^ esc" ^ ^ 

4+5 sin" 1 (|) 



92+ 


5tt 
75 sin" 1 ! 


i) 




75tt 




1 - 


2 






" 5V5 




924 


75 cos -1 


;§) 


V5- 


80^ 


77s) 



3tt 



0.785398 
0.707107 
0.833333 
0.883883 
0.577350 
0.816497 
0.746830 
0.863648 
0.979236 
0.886838 
0.879838 
0.549815 
0.685627 
0.821115 
0.903076 

0.504975 



\ 



(97) 



It might be of interest to address separability problems that appear to be computationally 
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intractable in the generalized Bloch representation by transforming them into the Bloore 
parameterization. 

We leave the reader with the intriguing questions, still left unanswered: do simple, exact 
brmulas exist for the Hilbert-Schmidt — and/or Bures, Kubo-Mori, Wigner-Yanase, . . . (cf. 
9j) — separability probabilities for the full 9-dimensional and 15-dimensional convex sets of 
real and complex 4x4 density matrices (and the 20-dimensional and 35-dimensional convex 
sets of real and complex 6x6 density matrices)? Further, are there helpful intuitions - 
above and beyond the direct implementation of our formulas ({TBI) and (ITT)) — that can aid 
in understanding the abundance of elegant, simple results occurring in this general research 
area? 
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